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Abstract 

We investigate the problem of reconstructing a fully anisotropic conductivity tensor 7 
from internal functionals of the form Vu • 7VU where u solves V • (7V11) = over a given 
bounded domain X with prescribed Dirichlet boundary condition. This work motivated by 
hybrid medical imaging methods covers the case n > 3, following the previously published 
case n = 2 [21]. Under knowledge of enough such functionals, and writing 7 = T7 (det7 = 
1) with t a positive scalar function, we show that all of 7 can be explicitely and locally 
reconstructed, with no loss of scales for r and loss of one derivative for the anisotropic 
structure 7. The reconstruction algorithms presented require rank maximality conditions 
that must be satisfied by the functionals or their corresponding solutions, and we discuss 
different possible ways of ensuring these conditions for C 1 '"-smooth tensors (0 < a < 1). 

1 Introduction 

Hybrid medical imaging methods aim to combine a high-resolution modality (such as acoustic 
waves or Magnetic Resonance Imaging) with a high-constrast one (e.g. Electrical Impedance 
Tomography, Optical Tomography, . . . ) in order to improve the result of the latter thanks to a 
physical coupling. In this context, the problem we consider is motivated by a coupling between 
an elliptic equation (modelling conductivity or stationary diffusion) and acoustic waves. Namely 
we consider the problem of reconstructing a fully anisotropic conductivity (or diffusion) tensor 
7 over a domain of interest X C W 1 from knowledge of a certain number of power density 
functionals of the form 1-L~[u](x) = Vu • yVu(x), where u solves the following partial differential 
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equation 



V • ( 7 Vu) = -J2 dih ij d jU ) = (X) 



u\ax = 9, 



(1) 



where the boundary condition g is prescribed. By polarization, we will see that mutual power 
densities of the form Vu • 7Vt> will also be considered, where both u and v solve (1) with 
different boundary conditions. The model above, when considered as a diffusion model for pho- 
tons in tissues, should be augmented with a term o a u accounting for absorption and will be 
addressed in future work. The availability of such functionals is justified by a coupling with 
acoustic waves, as it is described in the context of Ultrasoud Modulated- Electrical Impedance 
Tomography (UMEIT) or Optical Tomography (UMOT) in [2, 7, 15, 5] by considering acous- 
tic deformations, or in the context of Impedance-Acoustic Computerized Tomography in [13] 
by considering thermoelastic effects. In both cases, the acoustic waves come to the rescue of 
an otherwise very ill-posed problem (the classical Calderon's problem of recovering 7 from its 
Dirichlet-to-Neuman operator, see [10]), by providing internal functionals instead of boundary 
ones. In (1), we require 7 to have bounded components and to be uniformly elliptic as defined 
by the following condition 



from some k > 1. Borrowing notation from [3], we denote CX7) the smallest such constant n 
and define the set 



With these definitions, our problem may be formulated as follows 

Problem 1.1 (Inverse conductivity from power density functionals). For 7 in T,(X) or any 
subset of it, does the power density measurement operator % 7 uniquely characterize 7 ? If yes, 
how stably ? 

The problem just described has received fair attention in the past few years. The first inver- 
sion formula for Problem 1.1 was given in [11] in the isotropic, two-dimensional setting. There, 
a constructive algorithm as well as an optimal control approach for numerical reconstruction 
were presented. [16] then studied a linearized, isotropic version of Problem 1.1 in dimensions 
two and three with numerical implementation. 

Problem 1.1 has also been studied under constraints of limitations on the number of power 
densities available, the most restrictive case being the reconstruction of an isotropic tensor 
7 = al n in (1) from only one measurement H = o~\S/u\ 2 . In this case, a may be replaced in (1) 
by H/\\7u\ 2 , and this yields the following non-linear partial differential equation 



l^l 2 ^ 1 <7(xK^<K|e| 2 , xEX,£eM ; 



n 



(2) 



S(I):={7ei°°m, C( 7 )<oo}. 



(3) 
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Newton-based methods were proposed in [13] in order to successively reconstruct u and a, and 
the corresponding Cauchy problem was studied theoretically in [4]. 

In search for explicit reconstruction formulas using larger numbers of functionals, the authors 
first extended the reconstruction result from [11] to the three-dimensional, isotropic case in [5] 
with Bonnetier and Triki. This result was then generalized in [22] to re-dimensional, isotropic 
tensors with more general types of measurements of the form <t 2q |Vu| 2 with a not necessarily 
\. This covers the case a = 1 of Current Density Impedance Imaging [23, 24]. Finally, the same 
authors derived reconstruction formulas for the fully anisotropic two-dimensional problem and 
validated them numerically in [21]. 

In the last three papers presented, the explicit reconstruction algorithms were derived in the 
case where the power densities belong to W 1,00 (X), and assuming some qualitative properties 
satisfied by the solutions. In particular, the reconstruction algorithm for the isotropic case (or, 
equivalently, of a scalar factor multiplied by a known anisotropic tensor) strongly relies on the 
existence of re solutions of (1) whose gradients form a basis of K n at every point of the do- 
main. Under such assumptions, stability estimates were derived for the reconstruction schemes 
proposed, of Lipschitz type for the determinant of the conductivity tensor under knowledge of 
the anisotropic structure 7 := (det7)~~7, and of (less stable) Holder type for the anisotropic 
structure 7. Finally, it was shown for certain types of tensors 7 that the assumption of lin- 
ear independence made on the solutions could be guaranteed a priori by choosing appropriate 
boundary conditions, so that all the reconstruction procedures previously established could be 
properly implemented. 

Studying a linearized version of Problem 1.1 from the pseudo-differential calculus standpoint, 
the Lipschitz stability mentioned above was also pointed out in [17] in the isotropic case. There, 
the authors showed that from three power densities functionals, the linearized power density 
operator is an elliptic functional of an isotropic tensor a. They also studied in more detail the 
"stabilizing" nature of internal functionals of certain kinds that have arisen in all the hybrid 
medical imaging methods mentioned above. An extension of this result to the anisotropic case 
is presently investigated by the authors with Guo in [6]. 

The present work aims at unifying and extending the work done in [11, 5, 22, 21] by treating 
in full extent the anisotropic, re-dimensional case of Problem 1.1 for C 1,Q -smooth conductivities 
with < a < 1 (the Holder exponent is required by forward elliptic theory). The basis of this 
work also appears and will strongly rely on the first author's recent thesis [20]. 



2 Statement of the main results 

We decompose the conductivity tensor 7 into the product of a scalar factor r := (det7)» and a 
scaled anisotropic structure 7: 

7 := T7, r := (det7)«, det7 = l. (4) 
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Note that when 7 6 T,(X), r is uniformly bounded above and below by C("y) and C(~f)~ l , 
respectively. 

Under knowledge of enough power densities inside the domain, the reconstructibility of r 
and/or 7 are local questions, since under certain conditions described below, both quantities r 
and 7 can be explicitely and locally recovered in terms of power densities and their derivatives. 
We first describe these conditions and the corresponding recosntruction formulas in the next 
paragraph. 

Second, as the reconstruction algorithms presented above require local conditions, we will 
describe how to control these conditions from the domain's boundary, also tackling the question 
of global reconstruction. 

2.1 Local reconstruction algorithms 

Reconstruction of the scalar factor r knowing 7; We first consider the question of 
local reconstructibility of the scalar factor r under knowledge of the anisotropic structure 7. 
The main hypothesis here is that we use the mutual power densities Hij := Vuj • j'Vuj (for 
1 < i, j < n) of n solutions (u\, . . . ,u n ) of (1) whose gradients are linearly independent over a 
sub domain Q C X, a condition which we formulate as 

inf I det(V«i, . . . , Vn n )| > c > 0. (5) 

Under this assumption we are able to derive the following reconstruction formula: defining 
A = 72 to be the positive matrix squareroot of 7, and decomposing A into A = ^JtA with 
det^4 = 1, knowledge of 7 implies knowledge of A. Further defining Si := AVui for 1 < i < n, 
the data becomes := Si ■ Sj. Such vector fields satisfy the following PDE's 

d{A~ 1 Sif = dlogr A (A~ 1 Si) b and V • (ASi) = -V logr • ASi, 1 < i < n, (6) 

where the equality of two-forms expresses the fact that d(A~ 1 Sif = d 2 Ui = (exact forms are 
closed), and the scalar equality is deduced from the conductivity equation. Here the b exponent 
denotes the flat (or index-lowering) operator for the Euclidean metric. From these PDE's, one 
can derive the following formula, first established in [20, Lemma 4.3.1] as a generalization of 
earlier results in [5, 11, 21, 22]: 

Vlogr = -LHI - * (V(\H\*H jl ) ■ AS A Ar l S j = -Vlog \H\ + ~(VH jl ■ ASi)A~ x Sj, 
n V / n n 

(7) 

Equation (7) may thus be used to substitute Vlogr into the PDE's (6), and the resulting 
system becomes closed for the frame S = (Si, . . . , S n ). We then show that such a system may 
be rewritten as a first-order quasilinear system of the form 

VSi = Si(S, H, dH, A, dA), 1 < i < n, x e U, (8) 
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where Si is a Lipschitz functional of the components of the frame S. Here VSj denotes the 
total covariant derivative of the vector field Si, a tensor field of type (1, 1) that encodes all 
partial derivatives d p Sf. System (8) can thus be integrated over any curve to reconstruct the 
value of S from knowledge of S(xo) for fixed xq G £1. Once S is known throughout r can 
be reconstructed throughout £1 by integrating (7) in a similar fashion. The PDE's (7) and (8) 
are overdetermined and come with compatibility conditions which should hold as long as our 
measurements are in the range of the measurement operator. In such a case, this leads to a 
unique and stable reconstruction in the sense of the following proposition, first stated in [20, 
Prop. 4.3.6-4.3.7]: 

Proposition 2.1 (Local stability for logr). Consider two tensors 7 = tA 2 and 7' = r'A 2 in 

T,(X), where A and A' are known and with components in W 1,co (X). Let Q, C X such that 
the positivity (5) holds for two sets of conductivity solutions (u±, ■ ■ ■ ,u n ) and (u' l7 ■ ■ ■ ,u' n ) with 
respective conductivities 7 and 7' , call their corresponding data sets {Hij,H-A with components 
in W 1,co (X). Then the functions logr and logr' can be uniquely reconstructed with the following 
stability estimate 

|| logr - logr'|| w i,cx>(Q) < e + C (\\H - H'\\ w i,oc( X) + \\A- A'Wwi,™^)) , 

where the constant C does not depend on Q, and eq is the error committed at some xq € Vt. 

Such a stability statement shows that under the condition (5), the reconstruction of t\q is a 
well-posed problem in W 1,co (0,). Section 3.1 contains the proofs of equations (7) and (8). 

Reconstruction of the anisotropic structure 7, then of r: Here and below, we denote by 
M n (M) the space ofnxn matrices with its inner product structure (A, B) := AijBij = tr (AB T ). 
We now derive an approach to reconstruct the anisotropic structure 7 from additional measure- 
ments. We start from a basis of solutions (u±, . . . ,u n ) satisfying (5) over C X. Considering 
an additional conductivity solution v, we show that, although the solutions (tti, . . . ,u n ,v) are 
themselves unknown, the decomposition of Vv in the basis (Vui, . . . , Vit n ) is known from the 
power densities. Combining these equations with the PDE's satisfied by the solutions allows to 
derive linear orthogonality constraints on the product matrix AS (S here denotes the matrix 
with columns Si, ... , S n ). 

Thus, any additional solution v, by means of its power densities with the support basis, gives 
rise to a subspace V C M n (R) orthogonal to AS, moreover a basis of V is known from the data. 
The dimension of V is accurately controlled in [20, Prop. 4.3.8] and its maximal value is 

dim V < d M := 1 + n(n + 1) /2. (9) 

The matrix AS is arbitrary in M n (M) except for its determinant, known up to sign, thus AS 
requires n 2 — 1 independent constraints to be determined up to sign. This requires that we 
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consider enough additional solutions V\, . . . ,vi such that their corresponding spaces Vi, ■ ■ ■ ,Vi 
satisfy (i) dim^- =1 Vi = n 2 — 1, and (ii) fx^=i is spanned by a non-singular matrix (this 
condition should always hold true when measurements aren't noisy, as this orthogonal space is 
nothing but MAS). In mathematical terms, the proper condition to satisfy is as follows: for 
1 < i < I, let Mr^i, . . . , M(j) rfjv/ span Vi (they can be constructed from the data), and denote 

M := {M {i)j | 1 < i < I, l<j< d M }, #M = d M l, (10) 

rewritten more simply as Ai = {Mi | 1 < i < djtfl} below. Conditions (i) and (ii) mentioned 
above will hold if for i£S!, there exists an n 2 — 1-subfamily of At with nonzero hypervolume. 
With the notion of cross-product in Appendix A. 2, this condition may be written under the 
form 

inf (det(AA(/) J ff- 1 AA(/)))« > Cl > 0, (11) 

for some constant c±, where I(n 2 — 1, dj/l) denotes the set of increasing injections from [1, n 2 — 1] 
to [1, d][.jl] (i-e. / G X(n 2 — 1, dj^l) is of the form / = . . . , i n 2_i) with 1 < i± < ■ ■ ■ < i n 2_i < 
c?m0) an d where M(I) = A/'(Mj 1 , . . . , Mj 2 _ ) is the cross-product defined in Appendix A. 2. 
Under condition (11), we are able to reconstruct 7 and VlogT via formulas (40) and (42). This 
reconstruction is unique and stable in the sense of the proposition below. 

Proposition 2.2 (Local stability for 7 and logr). Consider two tensors 7 = T7 and 7' = t'j' 
in S(X). Let Q, C X where u±, . . . , u n , v±, . . . , v\ and u[, . . . , u' n , v[, . . . , v[ satisfy conditions (5) 
and (11). Then 7 and 7' are uniquely reconstructed from knowledge of the power densities of 
the above sets of solutions, and we have the following stability estimate 

||V(logr- logT')|U~ (f2) + ||7 - 7'IIl~(q) < C\\H - H'\\ w ^ m . (12) 

Remark 2.3. Although Proposition 2.1 required bounded derivatives on the anisotropic struc- 
tures 7, this is no longer the case here as the frame S is reconstructed algebraically instead of 
solving a differential system that involves derivatives of the anisotropic structure. This is in 
good agreement with the fact that the stability statement (12) is only stated in L°°-norm for 7. 

Remark 2.4. The scalar factor r is reconstructed with better stability than the anisotropic 
structure 7, for which there is locally a loss of one derivative. Although the reconstruction 
procedure presented was not yet proven optimal in terms of number of power densities involved, 
this loss of one derivative cannot be avoided and finds justification in the microlocal analysis of 
the linearized problem that will appear on future work [6]. 



6 



Local reconstructibility: In the light of the local reconstruction algorithms previously de- 
rived, a tensor 7 £ is locally reconstructible from power densities if for every x £ X, there 
exists a neighborhood f2 x 9 2 and n + / boundary conditions (g\, ... , g n , h\, . . . ,hi) such that 
the corresponding n first conductivity solutions satisfy condition (5) and the I remaining ones 
satisfy condition (11) (which then ensures via Proposition 2.2 that 7 is uniquely and stably 
reconstructible over Q x ). Based on the Runge approximation for elliptic equations [19], we then 
have the following generic result: 

Theorem 2.5 (Local reconstructibility of C 1,a tensors, a > 0). If ^ £ C l,a (X), then 7 is locally 
reconstructible from power densities. 

Remark 2.6. In a similar manner to [8], the proof of Theorem 2.5 is based on constructing 
solutions locally, that will fulfill conditions (5)-(ll), after which such solutions will be approxi- 
mated using the Runge approximation by solutions of (1) globally defined over X and controlled 
from the boundary. Although this result establishes local reconstructibility for a large class of 
tensors, the applicability remains limited insofar as the boundary conditions are not explicitely 
constructed. 

2.2 Global reconstructions 

The previous approach consisted in deriving explicit reconstruction algorithms under certain 
a priori conditions (linear independence or rank maximality) satisfied by a certain number of 
solutions of (1). These conditions may be checked directly on the power densities at hand. 
As the user only has control over boundary conditions in this problem, it is thus appropriate 
to define sets of admissible boundary conditions, for which the conditions mentioned above are 
satisfied globally. 

The first admissibility set is that of m-tuples of boundary conditions (m > n) such that, 
locally, n of the m solutions of (1) have linearly independent gradients. This ensures that the 
scalar factor r is uniquely and stably reconstructible throughout the domain. For 7 € T,(X), we 
call such an admissibility set G™ (m > n), subset of (H2(dX)) m , whose full definition is given 
in Def. 4.1. 

On to the global reconstruction of (7, r), Definition 4.2 constructs a second set of admissible 
boundary conditions. We first pick g £ Q™ for some m > n so that a basis of gradients 
of solutions may be available everywhere. Considering I > 1 additional solutions generated by 
boundary conditions h = (hi, . . . , hi), we say that h belongs to A™'\g) if the spaces (Vi, . . . , V/), 
generated by (v\, . . . , vi) as in the previous section, form everywhere a hyperplane of M n (M), so 
that 7 and r may be reconstructed with the stability of Proposition 2.2. 

While the construction of these sets is somewhat tedious, they allow us to define sufficient 
conditions for global reconstructibility of all or part of 7. In particular, they allow us to re- 
formulate a reconstructibility statement into a non-emptiness statement on sets of admissible 
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boundary conditions Q 7 or Ay, which are characterized by continuous functional of power den- 
sities. Namely, for a tensor 7 = 7-7, we have 

Theorem 2.7 (Global reconstructibility). 1. Under knowledge of a C 1 -smooth 7, i/ie func- 
tion t € VK 1,00 (X) is uniquely reconstructible if Q™ 7^ /or some m > n, with a stability 
estimate of the form 

|| logr - logr'Hu/i.oopQ < C (\\H - H'\\ w i,^ {x) + \\j - 7'|| w i,°°pf)) ■ (13) 

2. 7 is uniquely reconstructible if there exists m > n and I > 1 suc/i i/mf <7™ 7^ anc? 
«^7*' i (g) 7^ O where g € ura'i/i i/ie stability estimate 

Il7 - 7'IU°°(X) + II logr - logr'Htyi^^) < C\\H - H'\\ w i,oo( X) . (14) 

Combining this with the fact that the conditions of linear independence stated above can be 
formulated in terms of continuous functionals of power densities and their derivatives, we can 
deduce further useful facts about the sets £/ 7 and A-p all of which are established in [20, Sec. 
5.2.1], allowing us to draw the following conclusions (see Sec. 4.2): 

• The reconstruction algorithms presented above remain stable under C 2 -smooth perturba- 
tions of the boundary conditions. 

• Conductivity tensors that are close enough in C 1,a norm can be reconstructed from power 
densities emanating from the same boundary conditions. 

• The property of being reconstructible carries through push- forwards of conductivity tensors 
by diffeomorphisms, see in particular Proposition 4.3 below. 

With these properties in mind, global reconstructibility is thus established for conductivity 
tensors that are C 1,a -close to or push-fowarded from the cases below: 

1. If 7 = rl n for some scalar function r € H 2" +3+e with e > 0, then 7^ for n even and 
Q™ +1 7^ for n odd. The proof can be found in [22] and relies on the construction of 
complex geometrical optics solutions. 

2. For 7 = I n , straighforward computations (see Sec. 3.3 below) show that Id\gx £ 
and {xf — xf +1 }™~l\Qx £ A^' n ~ (Id\gx)- From this observation, one can cover the case 

of a constant tensor 70 by pushing forward the above solutions with the diffeomorphism 
_ i_ 

ty(x) = 7 2 x. 

Outline: The rest of the paper is organized as follows. Section 3 justifies the local reconstruc- 
tion algorithms. For the reconstruction of r, Section 3.1 provides proof of equations (7) and (8), 
Section 3.2 covers the reconstruction of all of 7, while Section 3.3 concentrates on proving The- 
orem 2.5. On to the question of global reconstructibility, Section 4 first studies the properties of 
the admissibility sets and A™' 1 before discussing what tensors are globally reconstructible. 
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3 Local reconstruction formulas 

3.1 Reconstruction of the scalar factor r 

Geometric setting and preliminaries: We equip X with the Euclidean metric g(U, V) = 
U -V = 5ijU l V J , where the Einstein summation convention is adopted. For x £ X, (e±, . . . , e n ) 
and (e , . . . ,e n ) denote the canonical bases of T X X and T*X, respectively. The flat operator 
coming from the Euclidean metric maps a vector field U = Wei to the one- form U = We 1 . 
We also denote by V the Euclidean connection, i.e. the Levi-Civita connection of the Euclidean 
metric, which in the canonical basis reads Vc/F = (U l di)V :> ej. 

Over a set f2 C X where (5) holds, we have the following decomposition formula, true for 
any vector field V over f2 

V = HP«(V ■ S p )S q , H*i = [H-%. (15) 

For any invertible symmetric matrix M, applying (15) to MV and multiplying by M _1 yields 
also the more general formula 

V = H pq (V ■ MS^M^Sq. (16) 

Proof of equation (7): The proof essentially relies on the study of the behavior of the dual 
coframe 1 of the frame (A~ 1 S\, • • • , ^l -1 ^). Let us denote = sgn(det(5i, . . . , S n )), constant 
throughout O by virtue of (5). Since S T S = H with S = [Si \ . . . \S n ], we have that 

det(S'i, ... ,S n ) = a VdetH = a detH^, H = {H P q}i< Pt q< n . (17) 
For 1 < j < n, let us define the vector field Xj by 

o-j := (-iy-\ (18) 

where the hat over an index indicates its omission. Xj is the unique vector field such that at 
every x G £1 and for every vector V £ T X Q, we have 

j 

Xj(x) ■ V = det^T 1 ^, . . . , A^Sj-u V , A~ l S j+u . . . ,^ _1 5 n ). 

In particular, we have that for any S , ^(M)-valued function M and any vector field V, 

MXj ■ V = Xj ■ MV = det^- 1 ^, . . . , MV, A _1 5 n ) 

= det Mdet((M" 1 A- 1 )5i, . . . , V, . . . , (M~ 1 A~ 1 )S n ), 

For (Ei, ■ ■ ■ , E n ) a frame, (yji, ■ ■ ■ ,cj n ) is called the dual coframe of E if u)i(Ej) = Sij for 1 < i, j < n. 



(A" 1 Si) b A • • • A {A~ l S.f A • • • A (A- 1 S n ) b 
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that is, we have that 



(MXjY 



<Tj det M -k 



A S n ) since we have, 



(M- L A- L Sif A--- /\(M- L A- L Sjf A--- /\(M-' L A- L S n y . (19) 

(Xi, ■ ■ ■ , X n ) is, up to some scalar factor, the dual basis to (A Si, • • • , 
for i ^ j 

Xj ■ A~ l Si = det(A- 1 5i, . . . , A~ l Si, A~ l S h A~ l S 

i j 

since the determinant contains twice the vector A~ 1 Si. Moreover, when i = j, we have 
Xj ■ A~ l Sj = det(A _1 5i, . . . , A _1 5 n ) = det A" 1 det (Si, ...,S n ) = a det [A^m 
where we used relation (17). Therefore we can use formula (16) with M = A" 1 to represent Xj 



as 



Xj = H kl (Xj ■ A- x S h )ASi = ao^W 1 dei^H^ASt. (20) 

k,l=l 1=1 

We now show that Xj is divergence-free, that is V • Xj = for 1 < j < n. Indeed, we write 



V ■ Xj = *d * X) = *d 



(A^S^ A • • • A (A^Sjf A • • • A {A^Snf 



0, 



since an exterior product of closed forms is always closed, thus we have 

V • Xi = 



] — u ) 1 < j < n. (21) 

Combining equations (20) together with (21), and using the identity V • (/V) = /V • V + V/ • V 
for / a function and V a vector field, we obtain 

= V ■ Xj = V • (H jl detiA^H^ASi) 

= det(A^m)VW l ■ ASi + W l V det(A^H^) • AS t + det{A~ l m)W l V ■ (ASi). 

The last term is zero since V • (ASi) = and the second term expresses the dotproducts of 
Vdet(yl~ 1 i72) with the frame A^S. Thus we use the representation formula (16) with M = A 
and divide by det(A~ 1 H^) to obtain 

Vlogdet^- 1 ^) = ^(Vlogdet^" 1 ^) • ASfiA^Sj = -(VH jl ■ ASfiA^Sj, 

which upon writing log det(A~ 1 H2) = — log det A + \ log det H yields 

V log det A = -V log det H + (VH jl ■ AS^A^Sj. 

We now plug in the rescaling A = t^A (so that det A = t% ), which implies A -1 = t~^A~ 1 , and 
notice that the terms involving r cancel out in the right-hand side of the last equation, thus (7) 
is proved. 
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Proof of Equation (8): We start by recalling the first equation of (6) 

d{A- l Sif = F b A {A~ l Si)\ l<i<n, F:=Vlogr, (22) 

where F is now considered a functional of (Si, . . . , S n ) and of the known power densities (this 
functional relation was obtained using the divergence equations in (6), which are of no further 
use here). The main tool to derive a first-order differential system for (Si, . . . , S n ) from (22) is 
Koszul's formula 

2(VuV) ■ W = Vu(V ■ W) +V V {U ■ W) -V W {U ■ V) - U ■ [V,W] - V ■ [U, W] + W ■ [U, V], 

(23) 

which expresses covariant derivatives in terms of dotproducts and Lie-Brackets of given vector 
fields. The dotproducts of Si, . . . , S n are known from the power densities, while the Lie Brackets 
[ASi, ASj] are known from (22), as the following calculation shows 



[ASj,j4Sj] — ASi 


■ VH jk 


-ASj 


• VH ik 


-d(A- l S k )\ASi,AS 3 


) 


= ASi 


VH jk 


-AS, 


■ VHik 


-F b A(2" 1 S fc ) b (lS i , 


AS,) 


= ASi 


■ Viljfe 


-AS, 


■ VH ik 


- H k ,F ■ ASi + H ki F 


■ ASj, 



(24) 

where we have used (22) and the characterization of the exterior derivative 

dU\V, W) = V V (U ■ W) - V W (U -V) — U-[V,W]. (25) 

Unless A = I n , the frames S and AS do not coincide, therefore one must modify formula (23) 
in order to obtain the promised system. Following [20], we first choose to represent the total 
covariant derivative of Si in the basis of tensors of type (1, 1) given by {Si ® (A~ 1 Sj) b }f j =l , in 
which the decomposition is explicitely given by 

VSi = H k m^(V XSq Si ■ S p ) S, ® (2~ 1 S fc ) b , (26) 

see [20, Lemma 4.3.4]. The subsequent work consists in analysing the term Si ■ S p , in 
particular removing all derivations on the Sj's by moving them onto either known data Hij or 
the anisotropic structure A. 

The first step is to establish the following "modified" Koszul formula 

2< ^AS q Si "> ' S P = ^ ' AS q Hi P + ^AS t H 1P ~ V ls p ^ ^ 
— [Si, Sp] • Sq — [Sq, Sp] • Si + [Sq, Si] • Sp, 

where we have defined [U, V] = — V^yU. Equation (27) is obtained in [20, Lemma 

4.3.2] by using the torsion-freeness and the compatibility of the connection with the metric. 



11 



The second step is to establish for any vector fields U, V the following relation 

[U, V] A = A~ l [AU, AV] - A A (U, V), (28) 

where A A is a vector- valued tensor of type (2,0), whose expression in local coordinates is ex- 
pressed as 

A A (U,V) = ^(U l V - V l U q )A- l [A h A q ], A x ■= A(;d t ), 

as established in [20, Lemma 4.3.3]. Therefore, the tensor A A encodes differential information 
about the anisotropic structure and is identically zero if A is constant. Plugging relation (28) 
into (27), and then using the known Lie brackets expression from (24), we arrive at 



LASj, • A Sq — [ASq, ASp] - A Si + [ylS^, ^ISj] • A Sp 

+ A A (Si, Sp) ■ Sq + A^(Sq, Sp) ■ Si — A A (Sq, Si) • Sp 
= ASq ■ VH ip + ASp ■ VH iq - ASi ■ VH pq + 2H pq F ■ AS* - 2H qi F ■ AS P 

+ A A {Si, Sp) • Sq + A^(Sq, Sp) • Si — A^(Sq, Si) • Sp. 

The last expression no longer differentiates the vector fields Si, which fulfills our goal. Plugging 
the last expression into (26) and simplifying expressions of the form (15), we arrive at the 
expression 

VSi = X - (S k ® V\ k + MJ lk ® (l~ 1 5 fc ) b + (ASi ■ VHJ k ) Sj ® {AT 1 S k ) [ 

+ (F ■ ASi)A~ l - AF (g> (A^Sif (29) 
+ ±H k miP(A A (Si, S p ) ■ S q + A A (S q , S p ) ■ Si - A A (S q , Si) ■ S p ) Sj ® (A^Skf, 

where we have defined the data vector fields 

U jk := (VH jp )HP k = -H jp (VHP k ), l<j,k< n. (30) 

The last thing to notice is that, with the expression of F given by (7), the right-hand side of 
(29) is a polynomial in the components of S of order at most five. Together with the a priori 
uniform estimate 

n 

El SAx) 1 2 < n max HsAx), 
xGX, Ki<n 

i=i — 

we deduce that VSj = Si(S,H,dH,A,dA), where Si is Lipschitz-continuous with respect to S 
so that the method of characteristic is a uniquely defined and stable. 
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3.2 Reconstruction of the anisotropic structure 7, then of r 

As explained in Section 2, we start with n solutions (u\,--- ,u n ) whose gradients satisfy the 
rank maximality condition (5) over some f2 C X. We will call (Vui, • • • , Vu n ) the support basis. 

Algebraic redundancies: Using the support basis above and formula (15), we have for any 
additional solution the following relation 

Vu (i ) = H^(Vu {i) ■ jVu p )Vu q = H™H [i)p Vu q . 

As a result, the power density of any two additional solutions un-\ and uu\ is computible via the 
formula 

H m) = Vu (4) • 7 Vu (j) = HPm (i)p H rl H [j)r Vu q • 7 Vu, = H? r H [l)p H {j)r , 

i.e. the mutual power density of any two additional solutions is algebraically computible from 
the mutual power densities of each of these solutions with the support basis. In other words, any 
additional solution uu\ adds at most n non-redundant dimensions of data, that is, the quantities 
H(j\ p for 1 < p < n. 

Algebraic equations for AS: Let us add an additional solution v = u n +i and consider the 
mutual power densities of these n + 1 solutions {-ffjj}i<jj< n +i. As explained in Appendix A.l, 
by linear dependence of n + 1 vectors in M n , one can find n + 1 functions 

Hi = (-l)* +n+1 det{ff p q}i< ?3 < n> l<g<n+l, q^i, 

IX = det{Hp q }i<p :q < n , 

such that 

n 

ViSi + [J,AVv = 0, (32) 

i=l 

where [i never vanishes over f2 by virtue of (5). In particular, the following expression is well- 
defined over f2 

n 

S v = -H^^fiiSi. (33) 

i=i 

We now apply the operators d(A~ 1 -)^ and V • (A-) to relation (33), and using the fact that 
d(A~ 1 Si) b = and V • (ASi) = 0, we arrive at the following relations 

z\ f\{A~ 1 S i ) = Q and Zi-ASi = 0, where Z { := V — . (34) 
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The first equation describes the vanishing of a two-form, which amounts to n(n — l)/2 scalar 
relations, obtained by applying this two-form to vector fields AS P , AS q for 1 < p < q < n: 

H iq Zi ■ AS P - H ip Zi ■ AS q = 0, 1 < p < q < n. 

Put in other terms and defining Z := [Z\ \ . . . \ Z n ] and S := [S± \ . . . \ S n ], these relations express 
the facts that the matrix Z T AS is traceless, and that the matrix HZ T AS is symmetric, which 
we may express as orthogonality statements of the form 

{AS, Z) = and {AS, ZHtt) = 0, f] £ A n (R). (35) 

In other words, the matrix AS is orthogonal to the following subspace of M n (R) 

V := {Z(\I n + HQ), (A,fi) G R x A n (R)}. (36) 

As established in [20, Proposition 4.3.8], we have that diml/ = 1 + r(n — (r + l)/2), where 
r = rankZ, with maximal value 1 + n(n — l)/2 when r £ {n — 1, n}. 

Reconstruction algorithm: Assume now that we have I > 1 additional solutions (v i, . . . , v{) 
generating spaces Vi, . . . , Vi of the form (36). Let {e p <g> e q — e q ® e p , 1 < p < q < n} be a 
basis for A n (R), then the space X)i=i ^» ^ s spanned by the following family 

= {Z,i, ZiH(e p ®e q - e q ®e p ) \ 1 <i <l, l<p<q<n}, #M = d M l, (37) 

with d,M defined in (9). Assuming that dim^' =1 Vi = n 2 — 1 throughout Q,, there is a n 2 — 1- 
family of A4 spanning it, from which we would like to reconstruct AS via a cross product 
formula, see Appendix A. 2. Now, for any n 2 — 1-subfamily (Mi, . . . , M n 2_i) of Ai, the cross 
product M{M±, ... , M n 2_i) is 

(i) either zero if (Mi, . . . , M n 2_i) is linearly dependent, 



(ii) equal to ± 



detA/ F (Afi,...,Af n a_i) 



det(AS) 

In the second case, we compute 



^45" otherwise. 



A/"(Mi, . . . , M n *_ x )H- x M*(Mi, M re2 _! 



detAT(Mi,...,M r 



2-1, 



2 



ASH~ 1 S T A T . 



det(AS) 

Using the fact that SH~ 1 S T = I n , A = A T , and {det{AS)) 2 = detH, we deduce the relation 
MH- l M T = {det{MH- l M T ))^^f, with M := M{M l ,...,M n 2^ l ). 
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This expression also covers the case (i), as the factor in front of 7 is zero if (Mi, . . . , M n 2_i) is 
linearly dependent. As we do not know a priori which subfamily of M is independent, we may 
sum over all possible cases. With the notation I(n 2 — 1,c£m0 introduced in Section 2, we sum 
the last reconstruction formula over all possible n 2 — 1-subfamilies of Ai 

Af(I)H~ 1 Af(I) T = 7, where T := Yl (^t(Af(I)H^Af(I) T ))^ . (38) 

l£l(n 2 -l,d M l) IeT(n 2 -l,d M l) 

J 7 is a sum of nonnegative terms which vanishes precisely when dim^ =1 Vj < n 2 — 1, so a 

way of formulating the fact that dim Yl\=i Vi = n 2 — 1 and that fX;!=i ^i) ^ s spanned by a 
nonsingular matrix is indeed 

inf F{x) > ci > 0. (39) 

When condition (39) is satisfied, 7 is uniformly and uniquely reconstructed over Q by the 
following formula 



7 



1 

J 



N{l)H~ l N{l) T , xen. 

IeT(n 2 -l,d M l) 

On to the reconstruction of r, we restart from (7) and rewrite it as 



n 



7Vlogr = - \H\~2 [V(\H\2W l ) ■ ASi) AS j: 



(40) 



(41) 



where 7 is obtained from (40). Again, we will use the cross-product expression to express V log r 
solely in terms of data. For / £ X(n 2 — 1, oImI), we have, 



detM(I) 



VdetH 



AS, 



where there is a sign indeterminacy. However, this indeterminacy disappears when considering 
expressions as in the right-hand side of (41): 



V(\H\2-W l ) -N{I)ei )M(I)e 



detM(I) 



V(\H\2H jl )-ASi)ASi 



VdetH 

(det(AA(/) J ff- 1 AA(/) T )^ (v(\H\^H jl ) ■ AS t ) ASj 



Summing over / G I(n 2 — 1,cImI), we arrive at 



Y (v(\H\^H jl )-M(I)eA AT(/) ej = T (V(\H\2-H jl ) ■ AS A ASj = F\H\^V \ogr, 



lei 
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which finally may be inverted as 



V logr 




x E CI. 



(42) 



This reconstruction formula guarantees a unique and stable reconstruction of r with no ambi- 
guity 

Proof of Proposition 2.2. The reconstruction of (7, r) is based on formulas (40) and (42). Putting 
definitions (31), (36) and (68) together, we see that the right-hand side of (40) is, at every point, 
a polynomial of power densities and their first-order derivatives. Since the denominator T in 
(40) is bounded away from zero, we clearly have a continuity statement of the form 



where the constant C degrades like c]~ with c\ the constant in (39). On to the reconstruction 
of log r, we can make the same observation as before judging by equation (42) and the fact that, 
since det7 = 1, the entries of 7 -1 are polynomials in the entries of 7. This leads to a stability 
statement of the form 



3.3 Proof of Theorem 2.5 

The proof of Theorem 2.5 uses the Runge approximation for elliptic equations, which by virtue 
of [19, Equivalence Theorem p. 442] is equivalent to the unique continuation property. The latter 
property holds for conductivity tensors with regularity no lower than Lipschitz [12]. The Runge 
approximation, as it is proved in [19, 8] for instance and adapted to our case here, states that 
if f2 CC X, then any function u £ H l {Q) satisfying V • (7V?/) = over Q. can be approximated 
arbitrarily well in the sense of L 2 (Q) by solutions of (1), provided that 7 is Lipschitz-continuous 
throughout X. In fact, we require a little more regularity here (7 £ C 1,a (X) with < a < 1) 
for forward elliptic estimates. 

Step 1. Local solutions with constant coefficients: Fix xq £ X and B-^ r = Bs r {xo) C X 



a ball of radius 3r (r tuned hereafter) centered at x$. Denote 70 := 7(^0) an d Aq := 7^. We 
first construct solutions to the problem with constant coefficients, whose power densities will 
satisfy conditions (5) and (11). Such solutions are given by 



l-l'\\ L ™ m <C\\H-H' 




W 1 '°°(X) 



□ 



u®(x) := x 1 — Xq, 1 < i < n, and for 1 < j < n — 1, 



u° n+j (x) := -(x - x ) ■ Qj(x - x ), Qj := A Q ^UjAq 



(43) 
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where we have defined Mj := ej ® e J — ej + i ® e 7 . These solutions satisfy V • (70 Vn) = 
throughout M n , and we trivially have 

det(Vu§, . . . , V<) = 1, x G M n , (44) 

so that condition (5) is satisfied throughout B% r . Moreover, condition (11) is also satisfied as 
direct calculations lead to Z% = Qi = AQ 1 MiA^ 1 for 1 < i < n — 1, and the matrix H° := 
{Vu^ 1 • 7oVu°}" :)=1 is nothing but 70. Thus the space of orthogonality is given by 

V = ^MQi + Qi 70 A n (R) = A^ 1 ( ^IRH; + UiA n (R) J ^Iq 1 . 

i=l \i=l / 

The last space between brackets can easily be seen to not depend on x and it spans the hyperplane 
of traceless matrices {In}" 1 ) so that V = {'yo} ± - In particular, condition (11) is satisfied for some 
constant c\ > independent of x. 



Step 2. Local solutions with varying coefficients: From solutions , we construct 

a second family of solutions {-uj} 2 ™^ 1 via the following equation 

V-( 7 V<) = (B 3r ), u r\ aB3r =ul l<i<2n-l, (45) 

thus the maximum principle implies that 

max ||<||z,°°(b 3t .) < 3r and max JK +j || L oo (B3r) < Cr 2 , (46) 

where the constant only depends on the constant of ellipticity C(7). The difference of both 
solutions satisfies, for 1 < i < In — 1, 

-V-( 7 V«-nO)) = V-((7-7o)V^) (5 3r ), «-^)| 9B3r =0, (47) 

where the right-hand side belongs to C°' a (B 3r ) with a uniform bound in < r < ro for some r$. 
Thus [14, Theorem 6.6] implies that 

IK - u °i\\c^{B 3r ) < C(\\u r i ~ «i lli-(Bsr) + ll F illci.«(B 3 f)) ^ C'hllc^CX), (48) 

where the first constant depends on n, C(^), W^Wc 1 ," an d -63,,. Interpolating between (46) and 
(48), we deduce the first important fact 

lin A , W u i ~ u iWc 2 (B 3r ) = °- ( 49 ) 

Remark 3.1 (Dependency of the constants on the domain). The constant in (48) depends on 
dB 3r , thus on r, however this dependency works in our favor when shrinking the domain. This 
can be seen by rescaling the problem x — > xq + rx',x' G -83(0) to keep the domain fixed, and 
studying the behavior of the constants w.r.t. the rescalings. 
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Step 3. Runge approximation (control from the boundary dX): Assume r has been 
fixed at this stage. By virtue of the Runge approximation property, for every e > and 
1 < i < 2n — 1, there exists gf G H^{dX) such that 

IK ~ u i\\L 2 (B 3r ) < e -> where u\ solves (1) with u\ \qx = g\- (50) 

Now applying [14, Theorem 8.24] using the fact that V • (7V(u? — u\)) = thoughout B^ r , we 
deduce that there exists [3 > such that 

IK - u i\\ce(B^) ^ C W u i ~ u i\\LHB 3r ) < Ce, (51) 

where the constant only depends on n, C(j) and r = dist (B2 r ,dB^ r ), in particular the same 
estimate holds with \\uf — u^\\i^,^B 2r ) on the left-hand side. Finally, combining (51) with [14, 
Corollary 6.3], we arrive at 

C C 

IK ~ u i \\c 2 (BT) - ^2 IK ~ U i \\L™(B 2r ) < ^2 £ 5 

where the constant only depends on a, n, C(7) and ||7||cv*(X)- Since r is fixed at this stage, we 
deduce that 

lin ?> iJSSP 1 IK - "illcafBr) = 0. (52) 

Completion of the argument: For any (lei, the following functionals are continuous 

C 1,Q! (0) x C 2 (ft) x C 2 (n) B (7, u, v) #(7, u, u) = Vu • 7V7; G W 1 ' 00 ^), 
[Ty i,oo (0)r (n+i)/2 9 {i T ii}l <^.< n _> det{fl-y}^ =1 G W 1 ' 00 ^), 

# = {flij}i<i<i<an-i -> F(H,VH) G L°°(n), 
where in the last case, T is defined in (38) with I = n — 1 and its the domain of definition is 

^l,oo^](2n-l)n thg condition mf det{^}" i=1 > 0. 

Step 1 established that det{H^}i<.i<j< n and T(H°,VH ) were bounded away from zero over 
B r . Due to the limits (49) and (52), there exists a small r > 0, then a small e > such that 
maxi<j<2n-i IK ~~ n i > llc 2 (B r (a;o)) * s so small that, by the continuity of the functionals mentioned 
above, det{Hfj}i<i<j< n and F(H £ , VH £ ) remain uniformly bounded from zero over B r , where 
we have denoted Hf- := Vuf ■ ^Vu e - for 1 < i,j < 2n — 1. Conditions (5) and (11) are thus 
satisfied over B r by the family {uf j 2 ™] -1 which is controlled by boundary conditions. The proof 
of Theorem 2.5 is complete. 
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4 Global questions 



4.1 Admissibility sets and their properties 

For compactness of notation, we denote by X(M, N) (M < N) the set of increasing injections 
from [1, M] to [1, N] (i.e. / G 1(M, N) has the form I = . . . ,i M ) with 1 < h < ■ ■ ■ < i M < 



The sets Q^: The first admissibility set is that of boundary conditions ensuring that the scalar 
factor t is uniquely and stably reconstructible. This requires the existence of, locally, n solutions 
with linearly independent gradients. Although one can easily choose m = n in two dimensions 
thanks to [1, Theorem 4], some counterexamples in higher dimensions [18, 9] show that one may 
need stricly more than n solutions in general, hence the definition below. 

Definition 4.1 (Admissibility set G™, m > n). Let 7 € S(X) be a given conductivity tensor. 

For m > n, an m-tuple g = (g\, .., <7 m ) £ {H^{dX)) m belongs to G™ if the following conditions 
are satisfied (denote Ui the solution of (1) with ui\qx = gi)' 

(i) The power densities Hij = Vui ■ 'yS/uj belong to W 1,co (X) for 1 < i, j < m. 

(ii) There exists a constant C g > such that 



Condition (i) above allows to construct a finite open cover of X in a generic manner, where 
to each open set can be associated a single basis of n solutions, see [20, Prop. 5.1.2]. This 
basis can then be used to reconstruct r throughout each 0^. Doing this for each {l^ and patching 
reconstructions appropriately allows to reconstruct r in a globally unique and stable fashion, as 
is summarized in [20, Theorem 5.1.4]. 

The sets A^: On to the global reconstruction of the anisotropy 7, we now define a second 
class of sets of boundary conditions, such that the solutions generated satisfy condition (11) 
throughout X. 

Let 7 such that G™ ^ for some m > n and pick g £ G™ with constant C g as in (53). 
By virtue of [20, Prop. 5.1.2], there exists an open cover made of balls O = {Qk)k=i °^ a 
constant C g > and an indexing function 1^ = (i(fc)i, • • • , i(k)n) £ ^( n > m ) f° r 1 < k < K 



N). 



inf V?\g](x) >C g , where OTg](x) := V det{H pg } MeI . 



(53) 




H(k) ■= {H pg , p,q e I(k)}, 



(54) 
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i.e. one may use {V?/j}j e / {fc) as a support basis over 0, k . Given an additional solution v a , we 
now construct over each £l k a basis for the space V based on the local support basis: 

V a \n k =M.Z a(k) + Z a{k) H {k) A n (R), where for 1 < j < n, 

■ill/ \ (55) 

Z a (k)ej ■= (-iy +n+1 V (det{H pq , p e 7 (fc)j g€ subs (I {k) ,j,a)} / detH (k) ) 

and where "subs (I(k)^(k)j^ a )" 1S obtained from 1/^ by replacing in^j by a. From a collection 
of I additional solutions, similarly to (37), we build over each £l k the family of matrices 

M\n k = {Z m , Z m H (k) (e p ®e 9 -e,®e p ) 1 < i < I, l<p<q<n}, (56) 

of cardinality cLmI with cLm defined in (9), so that we may rewrite it generically as 

M\n k = {M (k)i \l<i< d M l}- 

Definition 4.2 (Admissibility set .A™''(g) for g € G™)- For m > n, let us assume that g = 
(ffi >• >9m) S GJ?, and let (O = {Qk} k= i,I, Cg) an open cover, an indexing function and 
a constant associated to it. For I > 1, we say that I additional boundary conditions h = 
(hi,"' ,hi) G (Hz(dX)) belong to the set of admissibility A™' l (g) if there exists a constant 
Cg,h > such that the following condition holds 

in o f ^IsMlu^x) >C gjh , where (57) 

l<k<K x&l k 

^[g,h]|n fc := det(Af {k) (J)H- k ]Af {k) (J) T ^, (58) 

j£X(n 2 -l,d M l) 

M {k) (J) :=M(M {k)jl ,...,M (k)jn2 _ i ). 

With definitions 4.1 and 4.2 in mind, in the sense of the present derivations, we may say 
that a tensor 7 is globally reconstructible from power densities if G™ 7^ for some m > n and 

for g € 0™, .A™''(g) 7^ for I > 1 large enough. 

4.2 Properties of the admissibility sets 
Openness properties of Gj and A y : 

• For C 1,Q -smooth 7 an C 3 -smooth dX, the sets £? 7 and Ay are open for the topology of 
C 2 ' a (dX) boundary conditions ([20, Lemma 5.2.2]). 

• For C 1,Q -smooth 7 ([20, Lemma 5.2.3]). 
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Behavior of and with respect to push-forwards by diffeomorphisms: In the 

topic of inverse conductivity, diffeomorphisms are used in the anisotropic Calderon's problem 
to exhibit an obstruction to uniqueness. Here, these diffeomorphisms work in our favor in the 
sense that the property of being locally or globally reconstructible from power densities carries 
through push- forwards by diffeormorphisms. 

Let fy : X — > fy(X) be a H-^^-diffeomorphism where X has smooth boundary. Then for 
7 G T<(X), we define over fy(X) the push- forward of 7 by fy, denoted ty+'y, the tensor 

** 7 (y) := ({J^Dfy 7 Dfy T ) ofy~ 1 (y), y£ty(X), J* := det (59) 

As explained in [3], \l/*7 G Ti(fy(X)), and $ pushes foward a solution u of (1) to a function 
v = u o fy^ 1 satisfying the elliptic equation 

-V„ ■ (¥*7V„i;) = (*(X)), w| a( *(x))=^o*~ 1 , 

moreover ^ and ^lax induce isomorphisms of H l {X) and 5 onto and iJa (<9(^(X))), 

respectively. For our proofs based on pointwise estimates, we will add the further requirement 
that ^ satisfies a condition of the form 

C^ 1 < \ J^{x)\ < Ci},, x £ X, for some constant Cq, > 1. (60) 

We define the relation (7, -X") ~ {p/ 1 ,X') iff there exists : X — > X' a diffeomorphism onto X' 
satisfying (60), such that 7' = ^^7. It is clear that ~ is an equivalence relation. 
With these definitions in mind, our main observation is the following 

Proposition 4.3 (Prop. 5.2.4-5.2.5 in [20]). For 7 G S(X) and * : X -> #(X) a W 1 ' 2 - 
diffeomorphism satisfying (60), we /iawe /or any rn > n 

^ 7 = {g°r' : (61) 

Moreover, if g G £7™ /or some m>n, then we have 

A™:\(go^) = {ho^- 1 : hG^(g)}. (62) 

Remark 4.4. In other words, when a tensor 7 is reconstructible from power densities, then 
so is any tensor of the form fy+j with fy defined as above. Moreover, if (gi, . . . , g m , hi, . . . , hi) 
are boundary conditions on dX whose corresponding solutions allow to reconstruct 7 via the 
above explicit algorithms, then one may pick precisely (gi, . . . ,g m , hi, . . . , hi) o ^j^ 1 as boundary 
conditions on d(fy(X)) to reconstruct \l/*7. In particular, iffy fixes the boundary dX , then one 
may pick the same boundary conditions as 7 to reconstruct fy+j. 
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Proof of proposition 4-3: Let g € G™ for m > n. The corresponding solutions {u^^Li are 
being pushforwarded to functions Vi = Ui o over ^(X) whose power densities are denoted 
H[j = Vvj • [ x I/^7]Vi'j. For this proof, primed quantities will always indicate quantities referring 
to the push-forwarded problem. Using the chain rule and the definition of x I/*7, we have the 
transformation law 

H^x) = \Jy(x)\H' ij (V(x)), xeX. (63) 

Since the functions £>™[g] defined in (53) are homogeneous polynomials of power densities of 
degree re, we have the following relation 

V^[g](x) = \Mx)\ n V^ 7 [ g o^ l ](Ms(x)), xeX. (64) 

By virtue of condition (60), the left-hand side of (64) is uniformly bounded away from zero if 
and only if the right-hand side is as well, which concludes the proof of (61). 

On to the proof of (62), we first look at how things are being push- forwarded locally. As 
in the preliminaries before definition 4.2, an open over O = {Qk}k=i °f X yields an open cover 
{^(Qk)}k=i of ^(X) with the same indexing function /. This is because of the transformation 
law 

det^xU^, . . .,V x u in )(x) = Jy(x) detiVyV^,. . . , V y v in )(^(x)), x G X, 

which ensures that x Ui\i£i(k) ls a basis over fij. iff {^ y Vi} i£ j^ is a basis over \I>(Q fc ) with 
Vi = Ui o vj/" 1 . Using (63) and the chain rule, the matrices Z a n^ defined in (55) admit the 
transformation law 

Z a(k) = D* T Z' a(k) o * (Q k ). (65) 

In the definition (55) of the space V Q |ft fe , the scalar function |J$| appearing from the fact that 
= \J^,(x)\H'^(^(x)) may be absorbed by the space A n (M), so that we may write 

Thus the family M\n k (56), from the elements of which one construct cross-products, transforms 
as 

M\n k = DV T M'\ n n k) °* («*)■ 
Using formula (71), we deduce that for J £ X(re 2 — 1, #W4) and throughout 
M {k) {J)=M{M {k)ni ...,M {k)jn2 J 

= M{D^! T M[ k)n o % . . . , D* T M[ k) ^ o *) 
= ( J *TD*- 1 M(Ml k)jl ,...,M[ k)jn2 Jo*. 
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In particular, the function J^^fg, h]|n fe defined in (58) transforms according to the rule 



F?> l [EM\i 



I Jty 



2n-l~^- -r-m,l 



^[gor'.horV*. 



(66) 



Again, by virtue of (60), the left-hand side of (66) is bounded away from zero iff the right-hand 
side is bounded away from zero. Taking the minimum over 1 < k < K does not change this 
property, thus (62) is proved. □ 

A Linear algebra 

A.l Relations of linear dependence 

Lemma A.l. Let (Vi, . . . , V n+ i) be n + 1 vectors in M. n , and denote Hij = Vi-Vj for 1 < i, j < 
n + 1. Then the following linear dependence relation HiVi = holds with coefficients 



Hi 



det(Vi, . . . , V n ) ■ det(Vi, V n+1 , . . . , V n 



det{H pq | 1 < p < n, l<q<n+l,q^i}, 1 < i < 



(67) 



and n n+ i = det(Vi, . . . , V n ) 2 = det{Hjj-} 



ij fl<i,j<n- 



Proof. Define the /V s as m the statement of the function and let us show that Y17=i HiV, = 0- 
Consider the vector field defined by the following formal (n + 1) x (n + 1) determinant 



V = det 



/ Vl ■ V 1 ■■■ V\ ■ V n Vi • V n+l \ 

Vn-Vl ••• V n ■ V n V n ■ V n+1 
\ Vi ••• V n K+l / 



i.e. computed by expanding along the last row. Then we have 

n+1 



V = y^(~l)* ™ &<^{{H pq }l< p <n,l<q<n+l,qiti) Vi 
i=l 
n+1 

= ^(-l) i+ " +1 det(Vi, ...,V n ) det(Vi, • • • , V h . . . , V n+1 ) V 

i=l 

n 

det(Fi, . . . , V n ) ■ det(Vi, . . . , V n+1 , . . . , V n )V + det(Vl, . . . , V n ) 2 V n = ^ ^V u 



n+1 



i=l 



where moving V n+ i back to the i-th position in the i-th requires n — i sign flips. We now show 
that V = 0. For 1 < i < n, the dotproduct V ■ Vi becomes a determinant of a matrix whose 
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rows of indices i and n + 1 are equal, therefore V ■ Si = 0. Moreover, V ■ SVi+i is nothing but 
the determinant of the Gramian matrix of (Vi, . . . , V n +i), which is zero since n + 1 vectors are 
necessarily linearly dependent. Concluding, we have 

n+1 

V-V = ^2^V-Vi = 0, 

1=1 

thus V = 0, hence the lemma. □ 



A. 2 Generalization of the cross-product 

Let us consider a iV-dimensional inner product space (V, (, )) with a basis (ei, • • • , ejv)- Given 
a linearly independent family of N — 1 vectors (Vi,--- , Vat-i) in V, a (non-normalized) nor- 
mal to the hyperplane spanned by (Vi, • • • , Vjv-i) is given by computing the formal V-valued 
determinant 



TV(Vi, • • • ,Vat-i) 



1 



det(ei, • • • ,e N ) 



<*i,ei) 

(Vjv_i,ei) 
ei 



(Vi,e N ) 

(Vjv-i,ejv) 
e N 



(68) 



to be expanded along the last row. The function N can be easily seen to be N — 1-linear and 
alternating, and its definition does not depend on the choice of basis (ei, • • • , ejv)- Moreover, 
M satisfies the orthogonality property 

<A/-(Vi, • • • , Vn-^Vj) =0, 1 < j < N - 1, 

as such dotproducts take the form of determinants with identical j-th and iV-th rows. The 
first important property is that the squared norm of N represents the hypervolume spanned by 
Fi,...,Vat_i: 



(TV, TV) = det{(Vi, V^-)}i<ij<7v- 



(69) 



We now derive transformation rules when using linear transformations. For L : V — > V an 
automorphism, the following proposition relates TV(Vi, • • • , Vjy-i) with TV(LVi, • • • , LV^-i). 

Proposition A. 2. For L E Aut(V) and (Vi, • • • , V/v-i) a family of linearly independent vectors, 
the operator TV defined in (68) satisfies the transformation rule 



TV(L Vi , • • • ,LV]v_i) = (detL)L" T AA(Vi, 



(70) 
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Proof. Direct computations yield, picking a basis (ei, • • • , eyv 

{LVuei) ••• {LV u e N 

N{LV\, ■ ■ ■ ,LV N .i 



det(ei, • • • ,ejv) 



(LVSr-i,ei) 
ei 



(^Viv_i,ejv) 



detL 

det(L T ei, • • • , L T ejv) 



L 



-T 



(Vi,L T ei 



L T e 1 



{V u L T e N ) 



(V N -i,L T e N ) 



N 



where we recognize J\f(V\, • • • , Vjv-i) expressed in the basis (L T e±, • • • , L T e/v), hence the result. 

□ 

We are now interested in the case where V = Ai n (M.) with the inner product {Mi,M2) = 
tr (MiMg 1 ), and where the automorphism La denotes left-multiplication by a non-singular ma- 
trix A. First of all, it is straightforward to see that = L a t and L^ 1 = L^-i, where T and 
_1 on the right-hand side denote regular matrix transposition and inversion. 

With (e i, • • • , e n ) the canonical basis of W 1 , the family Eij — ej <S> for 1 < i, j < n is an 
orthonormal basis for M n (M) and we define the orientation on Ai n (M) by 

det(^n, • • • , E n i, • • • , E\ n , • • • , E nn ) = 1. 

Now, if we represent the vectors AEij in the above oriented basis, we see that 

AO 



det L A = det (AE U , ■■■ , AE nl , ■■■ , AE ln , ■■■ , AE n 



This brings us to the relation of interest: 



det 



'•• 

A 



Corollary A. 3. For (Mi, . . . , M n2 _i) G M n (R), A G G7 n (R) and M defined as in (68), we have 
the following transformation rule: 



J\f(AMi, . . . , AM n 2_x) = (det^) n ^- T AT(Mi,...,M n2 „ 1 ). 



(71) 
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